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COSMOLOGICAL EVOLUTION OF SUPERMASSIVE BLACK HOLES. I. MASS FUNCTION AT < Z < 2 
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ABSTRACT 

We present the mass function of supermassive black holes (SMBHs) over the redshift range z = 0-2, using 
the latest deep luminosity and mass functions of field galaxies to constrain the masses of their spheroids, which 
we relate to SMBH mass through the empirical correlation between SMBH and spheroid mass (the M, -Msph 
relation). In addition to luminosity fading of the stellar content of the spheroids, we carefully consider the vari- 
ation of the bulge-to-total luminosity ratio of the galaxy populations and the M,/Msph ratio, which, according 
to numerous recent studies, evolves rapidly with redshift. The SMBH mass functions derived from the galaxy 
luminosity and mass functions show very good agreement, both in shape and in normalization. The resultant 
SMBH mass function and integrated mass density for the local epoch (z « 0) match well those derived inde- 
pendently by other studies. Consistent with other evidence for cosmic downsizing, the upper end of the mass 
function remains roughly constant since z~2, while the space density of lower mass black holes undergoes 
strong evolution. We carefully assess the impact of various sources of uncertainties on our calculations. A 
companion paper uses the mass function derived in this work to determine the radiative efficiency of black hole 
accretion and constraints that can be imposed on the cosmological evolution of black hole spin. 
Subject headings: black hole physics - galaxies: evolution - quasars: general 



L INTRODUCTION 

Supermassive black holes (SMBHs) are ubiquitous in 
galaxies with a central stellar bulge (e.g., Magorrian et al. 
1998; Ho 1999; Kormendy & Kennicutt 2004), and, via the 
process of mass accretion (e.g., Salpeter 1964; Zel'dovich & 
Novikov 1964; Lynden-Bell 1969), power active galactic nu- 
clei (AGNs). The discovery of strong correlations between 
SMBH mass and the overall properties of the host spheroid, 
in particular the stellar luminosity, mass, and velocity disper- 
sion (Kormendy & Richstone 1995; Magorrian et al. 1998; 
Gebhardt et al. 2000; Ferrarese & Merritt 2000; Haring & Rix 
2004), has generated intense interest in the notion that black 
hole and host galaxy growth are closely connected. Study- 
ing the mass function of SMBHs and its evolution with red- 
shift is therefore of great significance not only for understand- 
ing their cosmological evolution of mass and spin (Shapiro 
2005; Wang et al. 2006, 2008, 2009), but also for addressing 
many issues in galaxy formation (e.g., Di Matteo et al. 2005; 
Shankar 2009). 

Generally speaking, there are two approaches to investigate 
the SMBH mass function (e.g., Tamura et al. 2006; Shankar 
2009). The first approach, which we refer to as theoretically 
based, stems from the pioneering work of Soltan (1982). 
The standard procedure, based on the argument that black 
hole growth is mainly driven by mass accretion (e.g.. Small 
& Blandford 1992; Marconi et al. 2004; Cao & Li 2008; 
Shankar et al. 2009b; Cao 2010), is to integrate the continuity 
equation to determine the number density of SMBHs and to 
couple it to the AGN luminosity function (LF). Such calcula- 
tions usually assume two parameters: the radiative efficiency 
for energy conversion of accretion and the Eddington ratio of 
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AGNs. As a result, comparing between the predicted SMBH 
mass function at z = with the locally observed quantity sets 
interesting constraints on the radiative efficiency. The radia- 
tive efficiency, averaged over redshift and black hole mass, 
is found to be ?7 « 0.1 (Yu & Tremaine 2002; Marconi et al. 
2004). 

The observational approach uses empirical scaling relation- 
ships between SMBH mass and host properties, in concert 
with the LF or stellar mass function (SMF) of the host galax- 
ies, to infer statistical information on the SMBH population 
(see Shankar 2009 for a review). Most previous works, be- 
cause of the availability of ample data, have focused on the 
local (z « 0) SMBH mass function (e.g., Graham et al. 2007; 
Vika et al. 2009; Shankar et al. 2009b). Few have treated 
the evolution of the SMBH mass function. Shankar et al. 
(2009a) predicted SMBH mass function at < z < 6 using 
the local velocity dispersion function of spheroids together 
with a model for their age distribution; however, their results 
are highly dependent on the assumed age distribution of the 
spheroids, which is poorly constrained. Out to intermediate 
redshifts (z « 1), Tamura et al. (2006) determined the SMBH 
mass function by employing a non-evolving relation between 
SMBH mass and spheroid luminosity (M, -Ljph relation) and 
by assuming that spheroids evolve passively. A number of re- 
cent studies, however, suggest that the black hole-host scaling 
relations evolve with redshift (McLure et al. 2006; Peng et al. 
2006a,b; Ho 2007; Bennert et al. 2010, 2011; Decarli et al. 
2010; Merloni et al. 2010). Although the exact magnitude of 
the effect is not yet well-defined, most studies find that the ra- 
tio of black hole mass to spheroid mass is higher in the past. It 
is essential to take into account these evolutionary effects. At 
the same time, outstanding progress has been made on deep 
multiwavelength surveys of normal galaxies and AGNs (e.g., 
Lawrence et al. 2007; Abazajian et al. 2009). This allows us to 
construct SMBH mass functions more reliably and to extend 
them to even higher redshifts than was possible before (e.g., 
Tamura et al. 2006; Vika et al. 2009, and references therein). 

In this work, we aim to derive the SMBH mass function us- 
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ing an up-to-date galaxy LF and SMF, with careful inclusion 
of various redshift-dependent effects. We extend and update 
the work of Tamura et al. (2006) out to z sa 2. It is interesting 
to emphasize that once the SMBH mass function is obtained, a 
combination of the above two approaches would immediately 
inform us of the redshift evolution of the radiative efficiency, 
from which inferences on black hole spins can be made. This 
is the subject of a companion paper 

This paper is organized as follows. In Section 2, we de- 
scribe the procedure for deriving the SMBH mass function. 
Section 3 presents the resultant mass functions and performs 
a comparison with previous works. We then investigate the 
uncertainties in our calculations and show how they impact 
upon the results (Section 4). We discuss the results of this 
work and the conclusions in Section 5. 

2. DERIVATION OF SMBH MASS FUNCTIONS 

We employ two methods to derive the SMBH mass func- 
tion: galaxy LFs, and, alternatively, galaxy SMFs. 

We use the recent /iT-band galaxy LF out to z w 4 obtained 
by Cirasuolo et al. (2010), who utilized the first data release 
of the UKIDSS Ultra Deep Survey. In addition to its depth 
and large survey area, the Ultra Deep Survey provides reli- 
able photometry in 16 bands. The K band is ideal for this 
application because it is not severely affected by dust absorp- 
tion and because it effectively traces the old stellar population 
and hence is a reliable indicator of stellar mass. 

Perez-Gonzalez et al. (2008) computed SMF for a sample 
of < z < 4 galaxies for which they could obtain rest-frame 
near-infrared photometry using Sp/fzer/Infrared Array Cam- 
era observations. They calculated stellar masses and photo- 
metric redshifts by fitting the observational data to ~2000 
reference templates for which there are reliable spectroscopic 
redshifts and well-covered spectral energy distributions from 
the ultraviolet to the mid-infrared bands. The stellar emis- 
sion models were computed using the code PEGASE (Fioc & 
Rocca-Volmerange 1997), adopting a Salpeter (1955) initial 
mass function (IMF) with stellar mass between 0.1 Mq and 
100 Mq. 

To derive the SMBH mass function from the galaxy LF or 
SMF, we need three ingredients: (1) the bulge-to-total lumi- 
nosity ratio (B/T) to transform the total galaxy luminosity to 
the spheroid luminosity, for which we assume that the light 
distribution follows the mass distribution in galaxies; (2) a 
scaling relation between SMBH mass and the mass of the 
bulge of the host; and, for the method involving the galaxy 
LF, (3) a prescription to describe the passive evolution of the 
spheroid luminosity to account for the evolution of the stellar 
population. Here we implicitly assume that all the galaxies 
under consideration possess a bulge component that hosts a 
black hole (Ho 2008), that the black hole mass scales with the 
bulge mass, and that the scaling relation evolves with redshift, 
as described below. 

Throughout the paper, we adopt a cosmological model with 
r^ni = 0.3, = 0.7, and Ha = 70km s"' Mpc"'. Unless other- 
wise specified, magnitudes are given in the AB system (Oke 
& Gunn 1983). According to Hewett et al. (2006), for the 
UKIDSS photometric system, the flux density for the zero 
point of the Vega-based magnitudes is f,^ = 6.31 x 10"^' erg 
s"' cm~^ Hz"'. The magnitude offset between the Vega- 
based system and the AB system is /STab = ^^veea + 1-9 (Hewett 
et al. 2006). We adopt a solar /iT-band luminosity of Lq.k = 
0.82 X 10^- ergs"' (see Table 2.1 in Binney & Merrifield 
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Fig. 1 . — Evolution of the bulge-to-total luminosity ratio, B/T. Data points 
show the mean B/T based on the type-dependent galaxy LFs of Zucca et 
al. (2006), and the en'or bars follow from the galaxy LFs. The solid line 
represents B/T = 0.6(1 +z)"'^ (see the text for details). 
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Fig. 2. — Passive evolution of the rest-frame A'-band absolute magnitude, 
calculated using PEGASE (solid lines) and GALAXEV (dotted Hues), with 
formation redshifts Zf = i, 4, and 10. For z < 1.5 the curves are well approx- 
imated by a straight line with a slope of ~ -0.7. 

1998), which corresponds to a solar /T-band absolute mag- 
nitude of Mq_k = 5.2 in the AB system. 

2.1. Bulge-to-total Luminosity Ratio 

B/T is an important morphological classification parame- 
ter that gives the fraction of the total luminosity contained in 
the bulge. Observations of local early-type galaxies show that 
they possess high values of B/T, while late-type galaxies pos- 
sess lower values of B/T (e.g., Simien & de Vaucouleurs 1986; 
Schechter & Dressier 1987). Many studies have performed 
bulge-to-disk decomposition of images of galaxies to quanti- 
tatively estimate their B/T (e.g., Kormendy 1977; Kent 1985; 
Aflen et al. 2006; Benson et al. 2007; McGee et al. 2008; 
Tasca & White 2011). Although its distribution is broad, 
B/T generally increases with galaxy luminosity (Schechter & 
Dressier 1987; Benson et al. 2007; Tasca & White 201 1). Be- 
low we will show how to determine the average value of B/T 
at z = 0, motivated by observations and the comparison be- 
tween our calculated SMBH mass function and the local ob- 
served one. 
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Fig. 3. — Mm -Lspi, relation in the K band. Lines sliow the M. — Ljph re- 
lation described by Equation (7) at j = (solid), 2=1 (dashed), and z = 2 
(dot-dashed). The values of the free parameters are listed in Table 1. Super- 
imposed for comparison are the observational data from Table 3 of Peng et 
al. (2006b; blue points) and Tables 2 and 3 of Bennert et al. (2010; black tri- 
angles and red squares). Note that the tabulated data from Peng et al. (2006b) 
and Bennert et al. (2010) were not corrected for luminosity evolution; the 
term accounting for the evolution of the spheroid luminosity in Equation (7) 
is essential to perform a direct comparison. We adopt l-con'ections appropri- 
ate for early-type galaxies, assuming K-V = —2.79 mag and K — R = —2.18 
mag (Fukugita et al. 1995; Gkardi et al. 2003). 



It is difficult to constrain B/T at high redshifts because of 
the challenges involved in obtaining reliable bulge-to-disk de- 
compositions for distant galaxies. However, we can use the 
fact that certain characteristic properties of galaxies statisti- 
cally correlate with B/T to estimate it indirectly at high red- 
shifts. Zucca et al. (2006), using template galaxy spectral en- 
ergy distributions, classified galaxies out to z w 1.5 into four 
spectral types — ranging from early-type to irregulars — and 
compiled morphological type-dependent galaxy LFs. Their 
four spectral types roughly correspond to the morphological 
types E/SO, Sa-Sb, Sc-Sd, and Irr. Based on the average prop- 
erties of nearby galaxies (e.g., Kent 1985; Weinzirl et al. 2009; 
Tasca & White 2011), we assign B/T w 0.7, 0.3, 0.1, and 0.0, 
respectively, to the above four type bins. Figure 1 plots the 
resultant variation of B/T with redshift, where the error bars 
follow from those of the galaxy LFs. The trend is physically 
reasonable: the mean B/T decreases systematically and mono- 
tonically with increasing redshift. The evolutionary trend can 
be parameterized as (see also Merloni et al. 2004) 



B/T = 0.6(1 -Fz)- 



(1) 



with 7 > 0, where the normalization is set to B/T = 0.6 at 
z = 0. As shown in Figure 1, 7 = 1 .3 gives a reasonably good 
description of the available data. 

This evolutionary trend of B/T can be understood within 
the framework of hierarchical galaxy formation (e.g., Croton 
2006; Khochfar & Silk 2006; Weinzirl et al. 2009), wherein 
spheroids are built up through two physical processes, merg- 
ers and secular evolution. Whereas major mergers scram- 
ble disks into (classical) bulges, multiple successive minor 
mergers can also transform a spiral galaxy into an ellipti- 
cal (Bournaud et al. 2007). Secular processes from global. 



non-axisymmetric structures (e.g., bars and spiral arms) in 
isolated galaxies further contribute to (pseudo-)bulge growth 
(Kormendy & Kennicutt 2004, and references therein). 

2.2. The M, — Msph and M,—Lsph Relation 

We parameterize the relations between black hole mass and 
spheroid mass and /T-band luminosity as 



log 



and 



M. 










M, 


= flLlog 









(2) 



(3) 



where {aM,bM,^M) and (a/,,/?/,, A/,) are parameters to be de- 
termined from the data, with Am and being the respec- 
tive intrinsic scatter of each relation. In our calculations, we 
adopt the M, - Msph relation derived by Haring & Rix (2004) 
and the M, -Lsph relation derived by Marconi & Hunt (2003), 
wherein (««, bM, Am) = ( 1 . 12 ± 0.06, -4. 12 ± 0. 10, 0.30) and 
(az„/'L,Az.) = (1.13 ±0.12, -4.11 ±0.07, 0.31). 

We emphasize that these two correlations pertain solely to 
nearby galaxies. Redshift-dependent effects should be consid- 
ered when applying them to higher redshifts. As mentioned 
above, mounting evidence suggests that the scaling relations 
between black hole mass and host galaxy evolve with redshift 
(McLure et al. 2006; Peng et al. 2006a,b; Ho 2007; Bennert 
et al. 2010, 2011; Decarii et al. 2010; Merloni et al. 2010). 
Broadly speaking, the ratio of the central SMBH mass and 
the stellar bulge mass, M,/Msph, tends to be higher at higher 
redshifts. If we express this redshift dependence as 

^•■oc(l+z)^ (4) 



M, 



sph 



with /3 > 0, the M, - Msph relation at any given redshift z can 
be rewritten 



M, Msph 

log — = AM log — — -I- /7m + P log( 1 + Z) . 

Mq Mq 



(5) 



Similarly, the M, - Lsph relation needs to be modified with the 
additional term /31og(l -l-z). 

Although most recent studies seem to be converging on the 
idea that black hole growth precedes bulge assembly at higher 
redshifts, there is no uniform consensus on the exact magni- 
tude of the effect (i.e., the value of /3). McLure et al. (2006) 
estimated the black hole-to-spheroid mass ratio in a sample 
of radio-loud AGNs, hosted by massive, early-type galaxies 
in the redshift range < z < 2, and found that /3 « 2. Peng 
et al. (2006b), analyzing the observed /?-band M, - Lsph rela- 
tion for a sample of 31 gravitationally lensed AGNs and 20 
unlensed AGNs at 1 < z < 4.5, concluded that, after account- 
ing for luminosity evolution, the M,/Msph at z > 1.7 is ^4 
times larger than the local value. The implied value of /3 is 
^ 1 .4. In a study of 89 broad-line AGNs between 1 < z < 2.2 
selected from the zCOSMOS survey, Merloni et al. (2010) 
inferred a somewhat weaker evolution for the ratio of black 
hole mass to host galaxy stellar mass, with /3 w 0.68; these 
authors did not explicitly separate the spheroid mass from 
the total galaxy mass, and it is possible that M,/Msph itself 
may evolve faster than indicated. A relatively strong degree 
of evolution in M,/Msph, consistent with (3 w 1.3-1.4, was 
reported by Bennert et al. (2010) for z ~ 0.4 and 0.6, by De- 
carii et al. (2010) for < z < 3, and by Bennert et al. (201 1) 
for 1 < z < 2. 
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TABLE 1 
Free parameters 
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Parameter 



-1 — I — I — I — I — I — I — I — I — I — r- 



1 galaxy LF 
• galaxy SMF 



Value 



Implication 



B/T 0.6 Ratio of bulge-to-total luminosity (B/T) at z = 

7 1.3 Power-law index for the red.shift evolution of B/T (Equation 1) 

/3 1.4 Power-law index for the redshift evolution of (Equation 4) 

™sph 

Q 0.7 Passive evolution of A'-band magnitude (Equation 6) 
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Fig. 5. — SMBH mass density (p.) as a function of redshift for M. > 
10^ A/0 derived from the galaxy LF (red squares) and the galaxy SMF (blue 
points). Redshift bins correspond to the bins chosen by Perez-Gonzalez et 
al. (2008). The black star symbol shows the locally observed value of p. = 
(4.9 ± 0.7) X 10^ Mq Mpc-'' (Vika et al. 2009). 



Fig. 6. — SMBH mass density (p.) as a function of redshift, normalized to 
z = 0, for Mm > 10^ Mq and for different combinations of /3— 7, derived from 
the galaxy LF (red squai'es) and the galaxy SMF (blue points). Redshift bins 
correspond to the bins chosen by Perez-Gonzalez et al. (2008). 



Based on the above summary, the weight of the cun^ent ev- 
idence suggests that M,/Msph evolves with redshift, and, for 
the parameterization given in Equation (4), that a reasonable 
choice for the redshift dependence of the evolution is /3 = 1 .4. 
We choose this as a fiducial value. Below we will show 
that the physical requirement that black hole mass density in- 
creases monotonically over time places additional constraints 
on (3. 

2.3. Passive Evolution of the Spheroid Luminosity 



When considering the M, -Lsph relation at high redshift — a 
necessary step for converting the galaxy LF into the SMBH 
mass function — a correction to the spheroid luminosity must 
be applied to account for the fading of the stellar population 
with time. (Note that our alternative procedure of using the 
galaxy SMF bypasses this complication, and thus serves as 
a useful consistency check for our LF-based method.) In 
practice, the correction depends on the star formation his- 
tory of the spheroid. Following common practice, we em- 
ploy a single-burst star formation model characterized by an 
e-folding time r and a formation redshift zj- Such a passive 
evolution scenario traces a single, maximally old stellar pop- 



COSMOLOGICAL EVOLUTION OF SMBHS . I. 5 



T u p 1 I 1 I 1 r-i-i-n-ni p 1 I 1 I 1 M 




M, /M^ M, /M^ M, /M^ 

Fig. 7. — Comparison of our SMBH mass function with tliat derived by Tamura et al. (2006) for different redsliift bins. Solid lines and dot-dashed lines are 
the mass functions derived from the galaxy LF and galaxy SMF, respectively, with hatched areas showing the errors. The shaded area denotes the mass functions 
derived by Tamura et al. (2006). 



ulation and gives a conservative estimate of the amount of 
luminosity evolution. 

We run the code PEGASE (version 2.0; Fioc & Rocca- 
Volmerange 1997) and GALAXEV (Bruzual & Chai'lot 2003) 
to compute the amount of luminosity evolution. Both sets of 
models use stellar evolutionary tracks from the Padova 1994 
library. We assume a Salpeter IMF, a lower mass limit of 
0.1 Mq, an upper mass limit of 100 M©, and solar metallic- 
ity. The e-folding time is set nominally to r = 1 Gyr, and 
the formation redshift is set to Zf = 3,4, and 10. The evolu- 
tion of the rest-frame /T-band luminosity is shown in Figure 
2. We find that the results are quite insensitive to zj. For com- 
pleteness, we also explore models assuming a ChalDrier (2003) 
IMF, as well as metallicities ranging from 0.4 Z© to 2.5 Z©, 
and T = 0.5-2 Gyr. The results are similar. In all cases, for 
z < 2, the luminosity evolution in the K band can be described 
as 

Mk{z) = Mk{z = 0)-Qz, (6) 

where Q w 0.7. Including this prescription for luminosity evo- 
lution, we can write the M, -Lsph relation at redshift z as 

log ^ = a, (log ^ - +b, + piog(l+z). (7) 

Figure 3 plots the M, - Lsph relation for the set of free pa- 
rameters tabulated in Table 1 . The observed data from Table 
3 of Peng et al. (2006b) and Tables 2 and 3 of Bennert et al. 
(2010) are superimposed for comparison. (Note that the orig- 
inal data were not corrected for luminosity evolution.) This 
comparison illustrates an important fact: a single M, - Lsph 
relation adequately describes the data from z = to z = 2 af- 
ter applying our prescriptions for the redshift dependence of 
M,/Msph and spheroid luminosity evolution. This verifies the 
robustness of our redshift-dependent corrections and reaffirms 



the conclusions of previous studies, such as those of Peng 
et al. (2006a,b), Bennert et al. (2010, 2011), and Decarii et 
al. (2010). 

Having introduced the above three components (for the red- 
shift evolution of B/T, M,/Msph, and bulge luminosity), we 
can relate the SMBH mass function to the galaxy LF, account- 
ing for intrinsic scatter in the M, - Lsph relation: 



$g(z,M.) = 2.5 I fl'logLsph.jf$(z,M^) 



dM, 



27rA, 



■ exp 



(logLsph./i:-logM.)2 



2Ai 



(8) 



where Mk is the total /iT-band absolute magnitude of the 
galaxy, Lsphx is the luminosity of the bulge, ^{z,Mk) is the 
A'-band galaxy LF given by Cirasuolo et al. (2010), and $g 
measures the mass function per unit volume per unit mass. 
Here the subscript "G" distinguishes the mass function of 
galaxies from that of AGNs. The factor dlogLsph.K/dM, = 
l/[ln(10)aLM,], according to Equation (3). The transforma- 
tion between Lsph,;^ and Mk reads 



MK=MQ_K-2.5\og 



1 i^sph.K 

B/T Lq,k 



(9) 



Expressing the galaxy LF as a Schechter (1976) function 
(Cirasuolo et al. 2010), 

$(2,M^) = 0.41n(10)$olO-"-^<"'"*'^'^exp[-10-"'4^'^] , (10) 

with AM = M^-M|, 

M*K{z)=M*Kiz = Q)-(^) , (11) 
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and 



5'o(z) = $o(z = 0)exp 



0.01 



(12) 



The free parameters are a = -1.07 ± 0.1, zm = 1-78 ± 0.15, 
;tM = 0.47 ± 0.2, = 1.70 ±0.09, /t^ = 1.47 ± 0.1, M*(z = 

0) = -22.26 (fixed), and $o(z = 0) = (3.5 ±0.4) x 10-^^Mpc"^ 
(see Table 1 in Cirasuolo et al.). 

Similarly, the SMBH mass function is related to the galaxy 
SMFas 



$g(z,M.) = 



dM, 

1 



-exp 



£/logMsph$(z,Mtot) 

' (logM,ph-logM.)2 



2A^ 



,(13) 



where Mtot is the total galaxy stellar mass, Msph = Mtot(B /T) 
is the mass contained in the bulge, and $(z,Mtot) is the 
galaxy SMF given by Perez-Gonzalez et al. (2008). The fac- 
tor c/ log Msph /(iM, = 1 /[ln(10)flMM,], according to Equation 
(2). The galaxy SMF, also expressed as a Schechter function 
(Perez-Gonzalez et al. 2008), is 



$(z,M) = ln(10)$*( ^ 



M 



(14) 



The parameters a, M*, and are given in Table 2 of Perez- 
Gonzalez et al. for 12 redshift bins from z = to z = 4. 

3. RESULTS 
3.1. SMBH Mass Functions 

Figure 4 shows the SMBH mass function, computed us- 
ing the galaxy LF and SMF, at z = 0, 1, and 2. The lowest 
redshift bin of Cirasuolo et al.'s LFs is at z = 0.3, and we ex- 
trapolate it to z = 0. We confirm that the extrapolation is con- 
sistent with the observed local /T-band galaxy LF (e.g., Dev- 
ereux et al. 2009, and references therein). The SMBH mass 
functions derived from the galaxy LF and SMF show strik- 
ingly good agreement, both in shape and in normalization. 
A notable discrepancy occurs at z = 2, where the SMF-based 
SMBH mass function systematically lies above that derived 
from the galaxy LF, forM, ^ (2-3) x lO^M©. As this corre- 
sponds roughly to the "knee" of the mass function, we expect 
this result to be quite sensitive to the values of and of 
the SMF. Curiously, inspection of the data in Perez-Gonzalez 
et al. (2008) reveals that M^, and show significant fluctu- 
ations precisely at the redshift bins z = 1.8 and z = 2.25; it is 
unclear whether the magnitude of these fluctuations is suffi- 
cient to account for deviations in the SMF-based SMBH mass 
function at z = 2. 

For comparison, we superpose the local (z ~ 0) SMBH 
mass functions derived by Marconi et al. (2004) and Vika 
et al. (2009). The agreement is excellent. We note that the 
mass function at the high-mass end hardly evolves from z ~ 2 
to z = 0, indicating that the most massive black holes have 
been largely in place since z ~ 2 and experience little growth 
since then (Marconi et al. 2004; McLure & Dunlop 2004). By 
contrast, the lower end of the mass function undergoes strong 
evolution. This trend is often referred to as cosmic "downsiz- 
ing." 

We calculate the mass density for SMBHs for M, > 
10^ Mq by integrating the mass functions and show the re- 
sults in Figure 5. Once again, the results derived separately 




Fig. 8. — SMBH mass function at ; = (top) and mass density as a function 
of z (bottom) for B/T = 0.4, 0.6, and 0.8. 

from the galaxy LF and SMF are quite consistent. The mass 
density increases rapidly from z = 2 until z ~ 0.5, below which 
it begins to saturate, converging toward the local (z ~ 0) mass 
density of p, = (4.9 ±0.7) x lO^MoMpc""' derived by Vika 
et al. (2009; black star). In light of the accretion paradigm 
for AGN activity and black hole growth, this is in line with 
the redshift evolution of AGN activity, which reaches a max- 
imum around z ~ 2-3 (e.g., Ueda et al. 2003; Hopkins et al. 
2007) and rapidly declines toward lower redshifts, such that 
by the present-day universe most black holes are highly sub- 
Eddington and nearly quiescent (Ho 2008, 2009). 

3.2. Constraints on the Free Parameters 

Table 1 lists the values of the free parameters we used to 
describe the redshift evolution of B/T, M,/Msph, and lumi- 
nosity of the stellar population. Their impact on the SMBH 
mass function can be inferred from Equations (7) and (9). In 
the Appendix, we explicitly show how these parameters affect 
the SMBH mass density. 

The preceding section explains how we adopt the values 
of key free parameters in our calculations, motivated by vari- 
ous lines of observations. Using the galaxy LFs for four bins 
of morphological types from Zucca et al. (2006), we obtain 
7 = 1.3 for the power-law redshift dependence of B/T and a 
normalization factor of B/T w 0.6 at z = 0. It is interesting 
to point out that an average B/T of 0.6 yields a black hole 
mass function in good agreement with the local observational 
data, as shown in Figure 4. Furthermore, that the SMBH mass 
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function derived from the galaxy LF matches well those ob- 
tained from the SMF testifies to the robustness of the spheroid 
luminosity evolution correction adopted in our calculations. 
Here we additionally show that the physical condition that the 
SMBH mass density always necessarily increases over cos- 
mic time imposes additional constraints on the evolution of 
the M,/Msph ratio. In the Appendix, we prove that the SMBH 



mass density derived from the galaxy LF is proportional to 

P,(Z) CX (1 +z)^"7«M0<-'/-'"''"/2-5""^ec/2.5^^(^)^ (^5) 

where is the slope of the M, - Lsph relation (see Equation 
(3)) and zm, ^m, and $o(z) are the parameters of the Cirasuolo 
et al. (2010) galaxy LF (see Equations (11) and (12)). We 
note that the parameters 7 and /3 are degenerate along lines 
of constant p-^ai. Since a/. « 1, for the purposes of this 
discussion let us set ai=\ for convenience. For a fixed 7 and 
a large value of /3, the first term in Equation (15) increases 
monotonically with redshift and dominates over the others; an 
increase of p, with redshift is physically prohibited. Figure 6 
illustrates p, derived for three choices of (i-^. While the case 
with /3-7 = -0.4 yields a rapid drop off of p, around z = 0, the 
case with l3-j = 0.4 leads to a negative increase. Both cases 
seem unlikely and conflict with the basic physical requirement 
that p, must increase with time and that most black holes are 
nearly quiescent by the present-day universe. Note that, for 
p, derived from the galaxy SMF, there is a sharp transition 
between the redshift bins z = 0.1 and z = 0.3; this may be due 
to the large uncertainties of the SMF at z = 0.1 (see Perez- 
Gonzalez et al. 2008). Inspection of Figure 6 suggests that, 
if 7 « 1.3, a reasonable value of /3 should be < 1.7. Future 
observations will empirically test this assertion. 

3.3. Comparison with Previous Results 

In Figure 7, we perform a comparison of our results with 
those of Tamura et al. (2006). We find good agreement for 
the massive end of the mass functions, for M, > 10^'^ M0, 
but the Tamura et al. results systematically fall below ours 
toward lower masses. We ascribe these discrepancies to three 
factors. 

(1) The survey depths of the galaxy LFs are very differ- 
ent. Tamura et al. employ the galaxy LFs from the COMBO- 
17 survey (Bell et al. 2004), whose depth in the B band is 
~ 25 mag (Vega), which corresponds to /iT w 21 mag (Vega). 
By comparison, the Cirasuolo et al. LFs reach K ^ 23 mag 
(Vega; Lawrence et al. 2007), allowing us to reach signifi- 
cantly further down the mass function, robustly sampling later 
type galaxies and lower mass black holes. 

(2) Tamura et al. only consider early-type galaxies, which 
occupy the luminous end of the galaxy LF and host exclu- 
sively massive black holes. We, on the other hand, have suffi- 
cient sensitivity to sample later type galaxies, which have pro- 
portionately lower mass bulges and hence black holes. Thus, 
it is not a surprise that our results agree with those of Tamura 
et al. at the high-mass end of the mass function but diverge to- 
ward lower masses. The excellent match between our results 
and those independently obtained by Marconi et al. (2004) 
and Vika et al. (2009) at z w (Figure 4) testifies to the ro- 
bustness of our assumptions. 

(3) Our calculations take into account the redshift- 
dependent evolution of B/T and M,/Msph, while Tamura et al. 
adopt constant values for these parameters. However, for a 
given total galaxy luminosity and our chosen prescription for 
passive evolution of the spheroid light,'* B/T and M,/Msph 
have opposite effects on the black hole masses. In other 
words, the parameters 7 and /3, which describe the evolution 
of B/T and M,/Msph, respectively, are degenerate. Because 
we use 7 = 1.3 and /3 = 1.4, in practice these two compo- 
nents nearly cancel each other out in terms of their effects 

Like us, Tamura et al. (2006) also compute passive evolution using the 
PEGASE code. 
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on black hole masses at high redshift (see Equations (7) and 
(9)). Therefore, inclusion of the redshift-dependent evolution 
of B/T and M,/Msph, in fact, induces only mild differences. 

4. SOURCES OF UNCERTAINTIES 

As the galaxy LP and SMF are well described by a 
Schechter function, the exponential decline beyond the knee 
implies that even minor changes in the high-luminosity or 
high-mass end can induce significant effects on the upper end 
of the SMBH function. Here we examine the impact of var- 
ious sources of uncertainties on our derivation of the SMBH 
mass function. For simplicity, we report only on results based 
on the galaxy LF; the results based on the galaxy SMF are 
similar 

4. 1 . Bulge-to-total Luminosity Ratio 

Figure 8 shows SMBH mass functions and mass densities 
for BAT = 0.4, 0.6, and 0.8 at z = 0. As to be expected because 
of the exponential drop off of the galaxy LF, the high-mass 
end of the black hole mass function is especially sensitive to 
the value of B/T. For a given total galaxy luminosity, larger 
B/T leads to larger black hole masses and hence to a larger 
number density and integrated mass density of SMBHs. The 



parameter 7, which measures the redshift evolution of B/T, is 
degenerate with respect to /3. Therefore, we treat ji-^ai as 
one independent parameter and have shown its influence on 
the SMBH mass density in Section 3.2 and in Figure 6. 

In our calculations, for the sake of simplicity, we neglect 
the dependence of B/T on the galaxy luminosity or mass, and 
just focus on its average dependence on redshift. Schechter 
& Dressier (1987) show that the bulge-to-disk ratio varies 
weakly with galaxy luminosity, roug hly asiL/L^f -^ (see also 
Figure 12 of Benson et al. 2007). For L « 0. IL* to lOL*, B/T 
increases from ~0.45 to 0.75, if we normalize B/T to 0.6 for 
L = Li,. This luminosity-dependent variation of B/T is small, 
and we neglect it. We further neglect any intrinsic dispersion 
in the average value of B/T at a given morphological type, as 
this dispersion is poorly quantified even for nearby galaxies 
(e.g., Simien & de Vaucouleurs 1986; de Jong 1996). 

4.2. The Local M, — Lsph Relation 

The upper panel of Figure 9 shows the fitting errors of the 
M, -Lsph relation derived by Marconi & Hunt (2003) and their 
corresponding influences on our results. We find that the re- 
ported uncertainties on the slope and normalization («£ and 
bi) of the M, - Lsph relation induce only minor effects on the 
overall SMBH mass function. Marconi & Hunt (2003) re- 
ported an intrinsic scatter of A/, w 0.3 for the M, - Lsph rela- 
tion. To assess its impact on our results, we calculate the mass 
function and mass density for three values of A/, (0.2, 0.3, and 
0.4), as shown in the middle and bottom panels of Figure 9. 
We find that A^ has a remarkably large impact on the shape 
of the upper end of the mass function. This is again due to 
the exponential decline of the galaxy LF. However, the effect 
of Ai on the mass density is relatively insignificant, changing 
the overall amplitude by ~ ±25%. 

Graham (2007) reanalyzed the work of Marconi & Hunt 
(2003), updating some of the data on black hole masses 
and galaxy magnitudes. He reports a somewhat shallower 
M, -Lsph relation, with best-fit parameters (aL,bL,Ai) = 
(0.93,-1.80,0.33). Figure 10 compares the effect that the two 
versions of the M, - Lsph relation have on our derived SMBH 
mass function, assuming an intrinsic scatter of Al = 0.3. The 
differences are minor over all masses. 

Two caveats should be kept in mind. First, recent stud- 
ies suggest that classical bulges and pseudo-bulges may obey 
somewhat different scaling relations with black hole mass 
(e.g., Greene et al. 2008; Jiang et al. 2011; Kormendy et al. 
201 1; Xiao et al. 2011). Another complication arises from the 
possibility that, at a given redshift, galaxies in different evo- 
lutionary phases may have different values of M,/Msph (e.g., 
Alexander et al. 2008; Lamastra et al. 2010). Unfortunately, 
the current data on the bulge properties of high-redshift galax- 
ies are far too crude to allow us to concretely address either of 
these issues. 

4.3. Choice of Passive Evolution 

We assume that the stellar population of the bulge formed 
in a single burst at a formation redshift of Z/=3, 4, and 10, 
thereafter fading passively as Mk{z) = Mk{z = 0) - Qz, with 
Q w 0.7. As we can see from Figure 2, the value of Q is quite 
insensitive to zy . In Figure 11, we examine the sensitivity of 
the SMBH mass density to the value of Q, choosing values of 
0.6, 0.7, and 0.8. Variations of this magnitude on Q have only 
a mild influence on p.. 

4.4. Summary of the Uncertainties 
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Among the uncertainties discussed above, the most impor- 
tant are those related with the local M, -Ljph relation. While 
uncertainties in the slope and zero point of the relation affect 
the SMBH mass function at the level of <; 0.2 dex over the 
entire mass range, the upper end of the mass function is espe- 
cially sensitive to the intrinsic scatter of the M, -Lsph relation. 
Changing the scatter at the level of ±0.1 dex affects the mass 
function by ±0.5 dex at M, « 2 x IO^'Mq. This is due to 
the exponential falloff of the galaxy LF, such that even a mi- 
nor change in luminosity leads to a significant variation in the 
galaxy number density, and hence in the SMBH number den- 
sity. In a similar manner, the upper end of the SMBH mass 
function is sensitive to the values of B/T. In our calculations, 
the redshift evolution of B/T, encapsulated in the parameter 7, 
is degenerate with the evolution of M,/Msph, which is param- 
eterized by /3. Fortunately, for reasonable choices of 7 and 
/3, which are well motivated by current observations, the red- 
shift evolution of B/T and M,/Msph largely cancel out. Plau- 
sible variations on the prescription for passive evolution of the 
spheroid luminosity also introduce relatively minor perturba- 
tions on our results. In summary, from inspection of the tests 
presented in Figures 6 and 8-11, it seems that the total un- 
certainties on the SMBH mass function are generally within 
-0.3 dex. 

Lastly, it is worth emphasizing that the stellar masses used 
to compute the galaxy SMF of Perez-Gonzalez et al. (2008) 
assume a Salpeter IMF. As has been pointed out by a num- 
ber of authors (e.g.. Bell et al. 2003; Bruzual & Chariot 2003; 
Pozzetti et al. 2007), the Salpeter IMF is too rich in low-mass 
stars and might overestimate the stellar masses. The Chabrier 
(2003) IMF may yield more reasonable stellar masses (e.g., di 
Serego Alighieri et al. 2005). Stellar masses derived with the 
Chabrier IMF are smaller than those derived with the Salpeter 
IMF by a factor of L7 (e.g., Pozzetti et al. 2007). Con- 
sequently, the characteristic mass M^, of the Perez-Gonzalez 
et al. galaxy SMF would be lower by this amount. In the 
power-law regime (M < M*), this is equivalent to a decrease 
of $i by a factor of — (1.7)'"""" (see Equation (14)). For a 
characteristic value of a « -1 .2 (Perez-Gonzalez et al. 2008), 
and thus the magnitude of SMBH mass function, de- 
creases slightly by ^ 0.05 dex in the power-law regime. Be- 
yond the knee (M > M^,), however, the exponential decline of 



the SMF significantly affects the magnitude of SMBH mass 
function. Nevertheless, the resultant SMBH mass density de- 
creases systematically only by a factor of — 1.7 (0.23 dex). 

However, the existing observations of the local SMBH mass 
function and mass density (e.g., Marconi et al. 2004; Vika et 
al. 2009) place constraints on the adopted IMF parameters. 
If the Chabrier IMF is adopted, a mismatch will be induced 
between the SMBH mass function and the existing observa- 
tions, and other parameters in the calculations will have to be 
modified accordingly (e.g., increasing B/T). Note that if the 
Chabrier IMF is universal and invariant with redshift, it will 
not affect the overall evolution of the SMBH mass function, 
and the present results remain unchanged. 

5. SUMMARY 

We derive the SMBH mass function from z = to z = 2 using 
the up-to-date, deep, wide-area /T-band galaxy LF of Cira- 
suolo et al. (2010), and, in a complementary manner, using 
the galaxy SMF of Perez-Gonzalez et al. (2008). In addition 
to extending much further down the mass function than pre- 
vious studies, enabling us to robustly sample later type galax- 
ies and lower mass black holes, our analysis carefully consid- 
ers redshift-dependent corrections to the average B/T of the 
galaxy populations, the M,/Msph ratio, and the stellar lumi- 
nosity of the bulge. We find excellent agreement between the 
SMBH mass functions derived from the galaxy LFs and the 
galaxy SMFs. Moreover, the resultant SMBH mass function 
and integrated mass density for the local epoch match well 
those derived independently by other studies. 

In a companion paper, we will use the SMBH mass function 
obtained in this work to determine the radiative efficiency of 
black hole accretion as a function of redshift, with the goal of 
constraining the cosmological evolution of black hole spin. 
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APPENDIX 

REDSHIFT DEPENDENCE OF SMBH MASS DENSITY 

In this appendix, we clarify how the SMBH mass density depends on certain free parameters adopted in our analysis. The 
Schechter (1976) function used to describe the galaxy LF has the general form 



$(z,L)flfL=$o(z) 



exp 



L\dL 



(Al) 



where ^o(z) is the normalization and L^, is the characteristic luminosity as a function of redshift. If, for the sake of simplicity, we 
neglect the intrinsic scatter of the M, -Lsph relation in Equation (3), the SMBH mass function ^q{z,M,) is related to the galaxy 
LF as 

^c{z,M,)dM, = '^>(z,L)dL. (A2) 



Integrating the SMBH mass function, we obtain the total black hole mass density 



$g(2,M.)M.JM. cx ^{z,L)M,dL 



According to the M, - Lsph relation (Equation 7), the black hole mass can be expressed as 

M, = lO'"--'"-^'''\l+zf[(B/T)L]'"-. 



(A3) 



(A4) 
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Hence, 

p,{z) = 10''^-"^Q^~/'~\1 +zf<^o(z)Uz){B/Tr'r{l-a + ai}, (A5) 

where is the slope of the M,-Lsph relation (see Equation (3)), $o(z) is given by Equation (12), and T{x) is the Gamma function. 
Since the characteristic luminosity from Equation (1 1) is 

Uz) = U{z = 0)10^^/^"''"/' ■^ (A6) 

we have 

P,{Z) OC (1 + .)'3-7«.10fcAM)'«/2.5-«,e-'/2.5^^(,)^ (^7) 

We note that, apart from being sensitive to the details of the galaxy LF, the redshift evolution of the SMBH mass density depends 
on the chosen prescription for the redshift evolution of B/T, the M, -Lsph relation, and luminosity variation of the stellar popula- 
tion. The parameters 7 and /3 are degenerate along lines of constant /S-'yaL- The physical requirement that p, always increases 
with time places additional constraints on the redshift dependence of these parameters (see also Hopkins et al. 2006). 
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